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CONVERSION  FACTORS,  METRIC  (SI)  U.  S.  CUSTOMARY 
UNITS  OF  MEASUREMENT 


Units  of  measurement  used  in  this  report  can  be  converted  as  follows 


Multiply 


BZ 


To  Obtain 


centimetres 

centimetres  per  millisecond 
grams  per  cubic  centimetre 
metres 

metres  per  millisecond 

millimetres 

megapascals 


0.3937007 

0.3937007 

62.42797 

3.280839 

3.280839 

0.03937007 

145.0377439 


inches 

inches  per  millisecond 
pounds  (mass)  per  cubic  foot 
feet 

feet  per  millisecond 
inches 

pounds  per  square  inch 


AN  ELASTIC-VISCOPLASTIC  CONSTITUTIVE  MODEL  FOR  EARTH  MATERIALS 


PART  I:  INTRODUCTION 

1.  The  mechanical  response  of  most  earth  materials  subjected  to 
high-intensity  dynamic  loadings,  such  as  those  produced  by  large  explo¬ 
sive  detonations  and  impact  of  high-velocity  projectiles,  differs  con¬ 
siderably  from  what  is  usually  observed  under  relatively  low-intensity 
static  loadings  (very  slow  rates  of  deformation) .  In  order  to  obtain 
realistic  solutions  for  dynamic  problems,  a  constitutive  model  should  be 
used  that  can  account  for  the  dependency  of  the  stress-strain-strength 
properties  of  earth  materials  on  the  various  rates  of  loading  or  defor¬ 
mation  being  applied.  In  addition,  the  constitutive  model  should 
account  for  other  pertinent  features  of  the  stress-strain  properties  of 
earth  materials  observed  under  both  dynamic  and  static  loading  condi¬ 
tions,  such  as  (a)  the  dependency  of  the  shearing  strength  of  the 
material  on  hydrostatic  stress,  (b)  shear-induced  volume  change,  and  (c) 
permanent  strain  during  a  load-unload  cycle  of  deformation  (for  both 
devlatoric  and  hydrostatic  loading  conditions) . 

2.  Incremental  elastic-plastic  constitutive  models  have  been  used 
successfully  to  simulate  the  stress-strain  properties  of  soil  (Baladi 
and  Rohani,  1979;  Baladi,  1977;  and  Sandler,  DiMaggio,  and  Baladi, 

1976) .  It  is  therefore  logical  to  adopt  a  physically  realistic  incre¬ 
mental  elastic-plastic  constitutive  model  for  earth  materials  and  intro¬ 
duce  rate  dependency  in  such  a  model.  Generally,  two  different  types  of 
rate-dependent  models  can  be  constructed  in  this  manner:  viscoelastic- 
plastic  models  in  which  both  the  elastic  and  the  plastic  responses  of 
the  material  are  rate-sensitive;  and  elastic-viscoplastic  models  in 
which  the  plastic  portion  of  the  model  is  rate-dependent  and  the  elastic 
portion  is  rate- independent.  As  pointed  out  by  Perzyna  (1966),  the 
viscoelastic-plastic  models  are  mathematically  very  complicated  and  are 
not  suitable  for  solving  practical  engineering  problems.  The  elastic- 
viscoplastic  models,  because  of  their  mathematical  simplicity 


(relative  to  viscoelastic-plastic  models)  and  their  similarities  with 
the  inviscid  theory  of  plasticity,  are  more  appropriate  for  practical 
engineering  application  (Perzyna,  1966;  Swift,  1975).  Also,  viscous 
effects  appear  to  be  more  evident  in  the  plastic  range  for  most  soils. 
Thus,  it  is  reasonable  to  adopt  an  elastic-viscoplastic  type  constitu¬ 
tive  relationship  to  model  the  rate-dependent  response  of  earth  mate¬ 
rials.  Application  of  this  type  of  model  for  describing  the  rate- 
dependent  response  of  rocks  is  discussed  by  DiMaggio  and  Sandler  (1971) 
3.  This  report  describes  the  fundamental  basis  of  elastic- 
viscoplastic  constitutive  relationships  and  the  development  of  a 
specific  model  of  this  type  for  earth  materials.  To  demonstrate  the 
application  of  the  model,  its  behavior  under  cylindrical  states  of 
stress  is  examined  and  correlated  with  experimental  data  for  a  clayey 
sand. 


PART  II:  FUNDAMENTAL  BASIS  OF 
ELASTIC-VISCOPLASTIC  CONSTITUTIVE  RELATIONSHIPS 


A.  The  basic  premise  of  elastic-plastic  constitutive  relation¬ 
ships  is  the  assumption  that  certain  materials  are  capable  of  undergoing 
small  plastic  (permanent)  as  well  as  elastic  (recoverable)  strains  at 
each  loading  increment.*  In  the  case  of  elastic-viscoplastic  materials, 
it  is  further  assumed  that  the  behavior  of  the  material  in  the  plastic 
region  is  rate-dependent.  Mathematically,  the  total  strain  rate  is 
assumed  to  be  the  sum  of  the  elastic  components  and  the  viscoplastic 
components ,  i . e . , 


de .  .  de  .  de  * 

ij  =  .  .11  +  _JJ 

dt  dt  dt 


where 


■  total  components  of  the  strain  rate  tensor 


.  E 

d£n 


dt  *  components  of  the  elastic  strain  rate  tensor 


— *  components  of  the  viscoplastic  strain  rate  tensor 

Elastic  Strain  Rate  Tensor 


5.  Within  the  elastic  range,  the  behavior  of  the  material  can  be 
described  by  an  elastic  constitutive  relation  of  the  type 


r  (  \ 

ijk£  dt 


*  In  this  paper,  compression  is  considered  positive. 
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material  response  function 


—  k£ 

— ^  ■  ■  components  of  stress  rate  tensor 
omr|  *  components  of  stress  tensor 

For  isotropic  compressible  elastic  materials ,  the  simplest  form  of 
Equation  2  is 


L.^l6  +L-—L J 

9K  dt  ij  2G  dt 


where 


J,  **  a  *  first  invariant  of  stress  tensor 

i  mm 


S. .  m  a. .  -  -r—  6 , .  =  stress  deviation  tensor 
ij  ij  3  ij 

K  -  elastic  bulk  modulus 

G  -  elastic  shear  modulus 


-f  “ 

lO  if 


*  Kronecker  delta  = 


6.  The  bulk  and  shear  moduli  can  be  expressed  as  functions  of  the 
Invariants  of  the  stress  tensor.  However,  in  order  not  to  generate 
energy  or  hysteresis  within  the  elastic  range,  the  elastic  response  must 
be  path-independent.  This  condition  can  be  met  if  and  only  if  the  bulk 
modulus  is  a  function  of  the  first  invariant  of  the  stress  tensor  and 
the  shear  modulus  is  a  function  of  the  second  invariant  of  the  stress 
deviation  tensor  (Sandler,  DiMaggio,  and  Baladi,  1976).  Thus, 


K  =  K(J.  ) ;  G  =  G(J1) 


(4) 


j;  ■  ?  su  sij 


the  second  invariant  of  the  stress  deviation 
tensor 


Viscoplastic  Strain  Rate  Tensor 

7.  The  behavior  of  the  material  in  the  plastic  range  is  assumed 
to  be  rate-sensitive.  As  indicated  by  Perzyna  (1966) ,  the  viscoplastic 
component  of  the  strain  rate  tensor  can  be  expressed  as  an  arbitrary 
function  of  the  "excess  stress"  above  the  initial  yield  condition  which 
is  called  the  static  yield  criterion.  The  static  yield  criterion  should 
satisfy  all  the  known  conditions  of  the  inviscid  theory  of  plasticity 
(Drucker,  1956).  In  general,  the  static  yield  surface  may  be  expressed 


<S(°ij*  K>  =  0 


For  isotropic  materials  the  static  yield  surface  may  be  expressed,  for 
example,  as 


<$s(Ji>  *0  =  0 


where  k  is  a  hardening  parameter  and  generally  is  a  function  of  the 
viscoplastic  (or  plastic)  deformation.  The  static  yield  surface 
(Equation  6)  may  expand  or  contract  as  k  increases  or  decreases, 
respectively  (Figure  1) . 

8.  Since  the  viscoplastic  strain  rate  is  an  arbitrary  function  of 
the  excess  stress  above  the  static  yield  criterion,  the  "dynamic  yield 
surface"  can  be  defined  as 


6d(Ji»  K»  = 


UJ1’  ^2>  K>  -  & 


and  the  flow  rule  for  work-hardening  and  rate-sensitive  plastic  mate¬ 
rials  may  be  written  as 


if  is>  0 

i£  (8i° 


(8) 


The  parameter  $/B  in  Equation  7  is  dimensionless  and  defines  the  rate 
of  expansion  of  the  yield  surface.  (B  is  a  scaling  factor  which  has 
the  dimension  of  stress;  it  is  introduced  in  Equation  7  to  make  ^ 
dimensionless.)  The  function  g(j{g/B)  in  Equation  8  is  a  dimensionless 
function  which  may  be  determined  from  the  results  of  dynamic  property 
tests  for  the  material  of  Interest  (Figure  2).  The  parameter  y  in 
Equation  8  is  a  viscosity  parameter  associated  with  the  viscoplastic 
response  of  the  material  and  has  the  dimension  of  (time)-'*'.  It  should 
be  pointed  out  that  for  very  slow  rates  of  deformation,  8  and  g(|{g/B) 
approach  zero.  Hence,  ^  and  become  identical  and  the  visco- 

plastic  flow  rule  (Equation  8)  reduces  to  its  inviscid  counterpart. 

9.  The  viscoplastic  stress-strain  relation  can  be  expressed  in 
terms  of  the  hydrostatic  and  deviatoric  components  of  strain.  Applying 
the  chain  rule  of  differentiation  to  the  right-hand  side  of  Equation  8 


The  deviatoric  component  of  the  viscoplastic  strain  rate  tensor 
de^/dt  can  be  written  as 


3  dt  ij 


Substitution  of  Equations  9  and  10  into  Equation  11  yields 


fiS.  /iis  !<s_ 
ic  2/JJ  3/jJ 


The  associated  flow  rule  is  satisfied  when  Equations  10  and  12  are  used 
with  Equation  7;  this  shows  that  the  vector  de^j/dt  in  the  (J^, 
space  is  always  normal  to  the  surface  »^»  *0  *  8  (Figure  3). 

Squaring  both  sides  of  Equation  9  gives 


SUBSEQUENT  DYNAMIC 
LOADING  SURFACE 


—LL 

dt 


INITIAL  YIELD 
'SURFACE  (<JS  =  0) 


Figure  3.  Dynamic  loading  surface  and  strain  rate  vector 


ll 


dt  dt 


Substitution  of  Equation  11  into  Equation  13  yields 


square  root  of  the  second  invariant  of  the 
viscoplastic  strain  rate  deviation  tensor 


Inverting  Equation  14  results  in 


I  fell)  +  2^ 

g-ll 1  3_\  dt  /  *Z\dt 

vji/  2Wn 


Comparison  of  Equation  7  with  Equation  15  yields  the  following  expres¬ 
sion  for  the  parameter  g/B 


Equation  16  Implicitly  represents  the  dependence  of  the  dynamic  yield 
surface  on  the  viscoplastic  strain  rate  tensor. 


Total  Strain  Rate  Tensor 


10.  The  total  strain  rate  tensor  can  be  obtained  by  combining 


Equations  3  and  9;  thus. 


PART  III:  ELASTIC-VISCOPLASTIC  MODEL  FOR  EARTH  MATERIALS 


11.  Before  the  development  of  the  elastlc-viscoplastic  model  is 
described,  it  may  be  useful  to  briefly  discuss  some  of  the  salient 
features  of  the  stress-strain  properties  of  soil  under  dynamic  loading. 
In  the  case  of  most  soils,  the  strength  of  the  material  Increases  with 
Increasing  rate  of  deformation.  Because  of  this  dependency  on  the  rate 
of  deformation,  the  overall  character  of  the  dynamic  stress-strain 
curves  often  differs  considerably  from  that  of  the  corresponding  curves 
obtained  under  static  loading  conditions  (Whitman,  1970) . 

12.  Figure  4  depicts  qualitatively  typical  stress-strain 
relations  under  axlsymmetric  triaxial  test  conditions  for  various  rates 
of  strain.  As  indicated  in  Figure  4,  the  soil  specimens  are  hydro¬ 
statically  consolidated  to  the  same  confining  pressure  (point  1,  Figure 
4).  The  samples  are  then  sheared  at  different  rates  of  strain  by 
increasing  the  vertical  stress  a z  while  the  radial  stress  o^  is  held 
constant.  The  important  behavior  to  be  observed  from  Figure  4  is  that, 
as  the  rate  of  strain  dez/dt  increases,  the  strength  of  the  material 
also  increases.  Therefore,  associated  with  each  strain  rate,  the 
material  possesses  a  unique  failure  envelope  which  may  be  referred  to  as 
"dynamic  failure  envelope."  The  implications  of  the  dynamic  failure 
envelope  can  best  be  realized  from  the  results  of  a  variable  strain  rate 
test. 

13.  It  is  possible  for  a  material  that  strain  hardens  under 
static  loading  conditions  to  exhibit  strain-softening  behavior  due  to 
strain  rate  effects  during  dynamic  loading.  For  example.  Figure  5 
depicts  the  hypothetical  result  of  such  a  test  superimposed  on  the 
corresponding  result  from  a  static  test.  Similar  to  Figure  4,  the  two 
stress-strain  relations  in  Figure  5  are  associated  with  axlsymmetric 
triaxial  tests  and  are  hydrostatically  consolidated  to  the  same  con¬ 
fining  pressure.  In  the  case  of  the  dynamic  test,  the  strain  rate 
during  the  initial  part  of  the  test  (point  1  to  point  2  in  Figure  5)  is 
increasing.  Beyond  point  2,  the  strain  rate  decreases.  During  the 
initial  part  of  the  test,  the  strength  of  the  material  continuously 


DYNAMIC  FAILURE  ENVELOPE 


Figure  4.  Typical  stress-strain  curves  for  different  strain  rates 
from  an  axlsymnetric  trlaxial  test 


Figure  5.  Dynamic  stress-strain  curve  for  a  variable  strain 
rate  from  an  axisymmetric  triaxial  test 


increases  because  of  the  increasing  strain  rate.  Beyond  point  2,  the 
strength  of  the  material  actually  decreases  because  of  the  decreasing 
strain  rate,  resulting  in  a  "falling"  or  softening  stress-strain  curve. 

14.  In  the  following  section,  the  mathematical  forms  of  the 
various  response  functions  are  developed  for  a  proposed  constitutive 
model  for  earth  materials. 


Elastic  Response  Functions 


15.  The  behavior  of  the  model  in  the  elastic  (recoverable)  range 
is  described  by  the  elastic  bulk  and  shear  moduli.  The  elastic  bulk 
modulus  is  assumed  to  be  a  function  of  the  first  Invariant  of  the  stress 
tensor  (Figure  6).  The  elastic  shear  modulus,  on  the  other  hand, 

is  assumed  to  be  a  function  of  the  second  invariant  of  the  stress 
deviation  tensor  (Figure  7),  i.e. , 


K 


[1  -  Kx 


exp(-K2J1)] 


(19) 


G 


exp(-G2/J^  )] 


(20) 


where 

K.  =  initial  elastic  bulk  modulus 

l 

and  K2  «■  material  constants 

G^  «  initial  elastic  shear  modulus 
G^  and  G 2  ■  material  constants 


The  constants  ,  and  K2  can  be  determined  from  the  charac¬ 

teristics  of  the  unloading  curve  from  an  isotropic  consolidation  test 
(Figure  8).  The  constants  G^  ,  G^  ,  and  G2  can  be  determined  from  the 


FIRST  INVARIANT  OF  THE  STRESS  TENSOR 


FIRST  INVARIANT  OF  THE 
STRESS  TENSOR,  J] 

Figure  6.  Elastic  bulk  modulus 
versus  first  Invariant  of  the 
stress  tensor 


Figure  8.  Proposed  relationship 
for  Isotropic  compression  test 


SQUARE  ROOT  OF  THE  SECOND 
INVARIANT  OF  THE  STRESS 
DEVIATION  TENSOR, 

Figure  7.  Elastic  shear  modulus 
versus  second  invariant  of  the 
stress  deviation  tensor 


Figure  9.  Proposed  relationship 
for  triaxial  shear  test 


a 


curves  from  triaxial  shear  tests  (Figure  9) . 

Viscoplastic  Behavior 

16.  For  the  plastic  behavior,  the  dynamic  loading  function  ^ 
(Equation  7)  is  assumed  to  be  isotropic  and  to  consist  of  two  parts 
(Figure  10):  a  rate-dependent  ultimate  failure  envelope  and  a  rate- 
dependent  strain-hardening  yield  surface.  The  failure  envelope  portion 
of  the  loading  function  is  assumed  to  be  of  the  Prager-Drucker  type  and 
is  denoted  by 


fd(Jl»  ^2*  e) 


f.«r  /5T>  -  8 


/jJ  -  oJ  -  k  -  B 


(2i; 


and  the  rate-dependent  strain-hardening  yield  surface  is  assumed  to  be 
elliptical  and  of  the  following  form 


Fd(Jl’  ^2*  K’  S) 


Fs(Jr  /jj,  <>  —  6 


V<v- 


(22: 


L)2/R2  +  J’  -  (k  -  L)/R  -  6 


B 


The  parameters  k  and  a  in  Equation  21  are  material  constants  rep¬ 
resenting  the  static  cohesive  and  frictional  strength  of  the  material, 
respectively.  The  parameter  R  (Equation  23  below)  in  Equation  22  is 
the  ratio  of  the  major  to  the  minor  axes  of  the  elliptic  yield  surface 
(Figure  10).  L  and  <  +  B  in  Figure  10  define  the  intersections  of 
each  elliptical  yield  surface  (Equation  22)  with  the  failure  envelope 
Fd^l*  ^2  *  ^  aru*  “*1  ax*s»  respectively.  The  hardening  parame' 

ter  k  is  generally  a  function  of  the  history  of  the  viscoplastic 
strain.  For  the  present  model,  R  and  k  are  chosen  to  be 


Strain-hardening  yield  surface 


b.  Rate-dependent  yield  surface 

Figure  10.  Proposed  yield  surface  for  the  elastic-viscoplastic  model 
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[l  +  R1exp(-R2L)] 
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(24) 


where  Rq  ,  R^  ,  R2  ,  D  ,  and  W  are  material  constants. 

17.  In  order  to  complete  the  specification  of  the  model  in  the 

viscoplastic  range,  the  function  g( <Jg/B)  (see  Equation  8)  must  be 

specified.  For  the  present  model,  g(^  /B)  is  chosen  to  be 

s 


on  the  failure  envelope 

■ 


on  the  hardening  surface 


(25) 


in  which  is  a  material  constant  controlling  the  rate  of  change  in 

the  value  of  the  function  g  .  Note  in  Figure  10  that  the  elliptic 
yield  surface  has  a  horizontal  tangent  at  its  intersection  with  the 
failure  envelope.  This  can  be  assured  by  the  following  relation: 


k  -  Rk 
L  “  1  +  oR 


(26) 


In  summary,  the  proposed  model  contains  16  material  constants.  Six  of 
the  constants  (K^  ,  ,  K2  ,  ,  and  G2)  are  related  to  the 

elastic  response  of  the  material,  and  the  remaining  10  constants  (B  , 
B^,a,k,Ro,R^,R2>D,W,  and  y)  describe  the  material's 
viscoplastic  response.  The  behavior  of  the  proposed  model  under  con¬ 
ventional  laboratory  test  conditions  (cylindrical  specimens)  is  examined 
in  the  following  section. 


PART  IV:  BEHAVIOR  OF  THE  MODEL  UNDER  A  CYLINDRICAL  STATE  OF  STRESS 


18.  The  rate-dependent  nature  of  the  proposed  model  can  be  demon¬ 
strated  by  examining  the  behavior  of  the  model  under  particular  labora¬ 
tory  test  boundary  conditions.  Since  most  mechanical  testing  of  soils 
for  engineering  purposes  is  performed  .ith  the  standard  triaxial  shear 
test  apparatus,  it  is  useful  to  investigate  the  behavior  of  the  model  in 
cylindrical  coordinates.  Adopting  the  z-axis  of  a  cylindrical  coordi¬ 
nate  system  (z  ,  r  ,  and  6)  as  the  axis  of  symmetry  of  the  sample,  the 
stress  tensor  and  the  strain  tensor  associated  with  this  configuration 
become 


and 


(27) 


(28) 


The  variables  P  =  J^/3  (mean  normal  stress),  (the  second  invariant 
of  the  stress  deviation  tensor),  and  3  (mean  volumetric  strain) 

associated  with  the  above  stress  and  strain  tensors  take  the  following 
forms 


J 


1 


a  +  2a 


Generally,  every  triaxlal  test  has  two  phases:  the  hydrostatic  phase 
and  the  shear  phase.  These  phases  are  discussed  below. 


Hydrostatic  Phase 

19.  During  the  hydrostatic  phase  of  the  triaxial  test,  all 
stresses  and  strains  are  equal.  Thus, 


a  =  a  =  J, /3  =  P 
z  r  1 


(32) 


e 
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(33) 


The  relationship  between  the  elastic  volumetric  strain  rate  and  the  rate 
of  the  first  invariant  of  stress  is  given  as  (Equation  3) 


de 


kk 


dt 


(34) 


in  which  the  elastic  bulk  modulus  K  is  given  in  Equation  19.  Sub¬ 
stituting  Equation  19  in  Equation  34  and  integrating  the  resulting 
expression  provides  the  following  relation  between  the  elastic  volu¬ 
metric  strain  and  the  first  invariant  of  the  stress  tensor: 


ili 

3K2K± 


iln 


exp(K2J1)  - 


(35) 


The  function  6/B  (Equation  16)  for  the  hydrostatic  phase  takes  the 
following  form: 


As  a  result  of  Equations  35  and  40,  the  total  volumetric  strain  for 
invi8cid  elastic-plastic  materials  takes  the  following  form: 
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kk 


1  -  K  |exp(K2J  )  -  K.l 

liqiq-  *"[ — nnq — J  +  *  - 


(41) 


As  expected,  the  volumetric  strain  in  Equation  41  is  independent  of  the 
time  rate  of  applied  loading. 


Shear  Phase 


21,  For  the  shear  phase,  a  constant  mean  normal  stress  test  will 
be  considered.  For  this  stress  path. 


dJ, 


dt 


=  0 


(42) 


de 


kk 


dt 


0  in  the  elastic  range 

vp 


(43) 


de 


kk 


dt 


in  the  viscoplastic  range 


de 


kk 


dt 


=  0 


(44) 


The  relationship  between  the  elastic  strain  deviation  rate  tensor  and 
the  rate  of  stress  deviation  tensor  is  given  as  (Equation  3) 


dt 


1_ 

2G 


(45! 


where  the  elastic  shear  modulus  G  is  given  in  Equation  20.  The  func¬ 
tion  S/B  (Equation  16)  takes  the  following  form: 


\(j1  -  l)  /r  + 


Equation  48  describes  the  coupling  between  the  volumetric  strain  and  the 
shear  stress  (J^  in  this  case)  during  a  constant  mean  normal  stress 
test.  For  a  specified  time  history  of  ,  Equation  48  can  be  inte¬ 
grated  to  determine  the  resulting  time  history  of  the  viscoplastic 
volumetric  strain.  From  Equations  12,  22,  and  25,  the  following  expres¬ 
sion  is  obtained  for  the  viscoplastic  strain  deviation  rate  tensor: 


L)/R  +  J 


yj1  -  L)  /R  + 


Equations  45  through  49  provide  a  complete  specification  for  the  behav¬ 
ior  of  the  material  for  constant  mean  normal  stress  triaxial  tests. 

These  equations,  however,  must  be  integrated  numerically  in  order  to 
relate  stresses  to  strains  during  dynamic  loading.  A  computer  program 
has  been  developed  for  numerical  integration  of  the  governing  equations 
of  the  proposed  model  for  general  three-dimensional  states  of  stress. 

For  the  sake  of  brevity,  the  numerical  Implementation  of  the  model  is 
not  described  herein.  In  the  following  section,  the  behavior  of  the 
model  is  correlated  with  test  data  for  a  clayey  sand  using  this  computer 
program. 


PART  V:  CORRELATION  OF  TEST  DATA  WITH  MODEL  BEHAVIOR 


22.  In  this  section,  the  behavior  of  the  model  under  states  of 

uniaxial  strain  and  triaxial  shear  are  correlated  with  available  test 

data  for  a  clayey  sand  classified  as  SC  according  to  the  Unified  Soil 

Classification  System.  Static  data  consisted  of  (a)  load-unload  axial 

stress-axial  strain  relations  (o  versus  e  )  for  uniaxial  strain,  and 

z  z 

the  corresponding  stress  paths  expressed  in  terms  of  principal  stress 
difference  versus  mean  normal  stress  (a  -  a  versus  P)  ,  and  (b)  two 

Z  IT 

triaxial  shear  test  load-unload  stress-strain  relations  (for  two 
different  confining  stresses)  presented  in  terms  of  principal  stress 
difference  versus  principal  strain  difference  ( a -  versus 
e  -  e  )  ,  and  a  static  failure  envelope  based  on  these  two  tests 

Z  IT 

(Ehrgott,  1978).  The  available  dynamic  data  for  this  material  consisted 
of  several  stress-strain  curves  from  dynamic  uniaxial  strain  tests 
(Jackson,  Ehrgott,  and  Rohani;  1980). 

23.  The  first  step  in  correlating  the  behavior  of  the  model  with 

test  data  is  to  simulate  the  static  properties  of  the  material  and  to 
determine  the  numerical  values  of  the  material  constants  , 

K2  ,  G±  ,  G^  ,  ,  a  ,  k  ,  Rq  ,  R^  ,  R2  ,  D  ,  and  W  .  The 

response  of  the  material  under  static  loading  is  governed  by  these 
constants  (i.e.,  the  case  of  inviscid  plasticity  when  rate  dependency  is 
neglected).  Figure  11  portrays  the  static  test  data  and  the  correspond¬ 
ing  model  behavior  for  both  the  triaxial  shear  and  the  uniaxial  strain 
test  conditions.  It  can  be  noted  from  Figure  11  that  the  model  has 
reasonably  simulated  the  static  response  of  the  soil  both  qualitatively 
(triaxial  shear  response)  and  quantitatively  (uniaxial  strain  response). 
The  next  step  is  to  simulate  the  available  dynamic  stress-strain 
properties  of  the  material,  using  the  numerical  values  of  the  13  con¬ 
stants  above  and  numerical  values  determined  for  the  remaining  three 
constants  B  ,  ,  and  y  .  As  indicated  before,  the  available 

dynamic  data  for  this  material  are  limited  to  several  uniaxial  strain 
stress-strain  relations.  The  dynamic  data  were  obtained  for  loading 


VERTICAL  STRESS,  a  ,  MPa  STRESS  DIFFERENCE 


•T 


VERTICAL  STRAIN,  ez>  PERCENT  MEAN  NORMAL  STRESS,  p,  MPa 

b.  Uniaxial  strain  response 

Figure  11.  Static  response  of  the  material  in  triaxial  shear  and  uniaxial 
strain  tests;  model  behavior  versus  test  results 


rise  times  on  the  order  of  a  few  tenths  of  a  millisecond  (Jackson, 
Ehrgott,  and  Rohani;  1980).  Two  of  the  dynamic  stress-strain  curves 
were  used  for  the  purpose  of  model  fitting.  Figure  12a  depicts  the 
stress-time  histories  for  the  two  tests.  It  should  be  noted  that  the 
entire  load-unload  cycles  for  these  tests  were  completed  in  slightly 
less  than  1  msec.  These  stress-time  histories  were  used  as  input  for 
driving  the  model.  The  experimental  stress-strain  curves  and  the 
corresponding  model  behavior  for  the  two  dynamic  tests  are  shown  in 
Figure  12b.  The  agreement  between  the  dynamic  test  data  and  the  model 
behavior  is  very  good,  both  quantitatively  and  qualitatively,  especially 
during  the  early  part  of  the  test.  It  is  of  interest  to  compare  these 
dynamic  stress-strain  relations  with  the  corresponding  static  result  in 
Figure  lib.  In  the  case  of  the  static  test,  the  stress-strain  curve  is 
concave  to  the  stress  axis  (a  "stiffening"  behavior),  whereas  the 
dynamic  curves  are  concave  to  the  strain  axis  (a  "yielding"  behavior). 
The  proposed  constitutive  model  predicts  this  dramatic  change  in  the 
overall  character  of  the  stress-strain  curves  remarkably  well. 

24.  The  numerical  values  of  the  16  material  constants  for  the 
clayey  sand  are  summarized  in  Table  1.  To  demonstrate  the  application 
of  the  proposed  model  further,  the  model  was  driven  with  the  strain-time 
history  shown  in  Figure  13a  under  triaxial  shear  test  conditions  (o^  •» 
2.07  MPa)  using  the  material  constants  given  in  Table  1.  The  resulting 
dynamic  stress-strain  curve  is  shown  in  Figure  13b.  Also  shown  in 
Figure  13b  is  the  corresponding  static  stress-strain  curve  (from  Figure 
11a).  The  effect  of  the  strain  rate  on  the  stress-strain  response  of 
the  material  is  clearly  demonstrated  in  Figure  13b.  As  the  strain  rate 
increases  during  the  early  part  of  the  simulated  dynamic  test,  the 
stiffness  and  strength  of  the  material  also  increases  (relative  to  the 
static  stiffness  and  strength).  During  the  latter  part  of  the  test 
where  strain  rate  decreases  with  time,  the  dynamic  curve  actually  falls 
(apparent  strain  softening)  and  eventually  coincides  with  the  static 
stress-strain  curve  (at  late  times). 


! 

Table  1 

Summary  of  the  Material  Constants  and 
Numerical  Values  for  Clayey  Sand 


Material 

Constant 


Unit  of 
Measure 


Numerical 

Value 


VERTICAL  STRESS,  o_,  MPa  VERTICAL  PRESSURE,  MPa 


VERTICAL  STRAIN,  e_,  PERCENT 


a.  Input  strain-time  history 


b.  Resulting  dynamic  stress-strain  relation 


Figure  13.  Predicted  dynamic  stress-strain  relation  for  triaxial 
shear  test  condition  (o  =  2.07  MPa) 


PART  VI:  SUMMARY  AND  CONCLUSIONS 


25.  An  elastic-viscoplastic  constitutive  model  has  been  developed 
for  earth  materials  and  has  been  partially  validated  via  comparison  with 
both  static  and  limited  dynamic  stress-strain  data  for  a  clayey  sand. 

The  model  is  capable  of  simulating  many  pertinent  features  of  the 
stress-strain-strength  properties  of  earth  materials  such  as  dependency 
of  the  shearing  strength  of  the  material  on  hydrostatic  stress  and  rate 
of  deformation,  shear-induced  volume  change,  and  permanent  deformation 
under  hydrostatic  and  deviator ic  cyclic  loadings.  In  its  present  form, 
the  model  contains  16  material  constants  that  can  be  readily  determined 
from  the  results  of  static  and  dynamic  triaxial  shear  and  uniaxial 
strain  tests. 

26.  The  model  has  been  translated  into  a  numerical  algorithm  for 
implementation  into  finite-difference  or  finite-element  computer  codes. 
The  numerical  algorithm  is  very  versatile  in  that  it  embodies  all 
classes  of  elastic-plastic  constitutive  models.  Test  data  for  several 
rates  of  deformation  and  test  boundary  conditions  (other  than  those  used 
to  fit  the  model)  are  needed  to  further  validate  the  accuracy  and  deter¬ 
mine  the  range  of  application  of  the  model. 
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